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We analyze extant data of Al 2+ , Al + and other low ions with the aim of studying the ioniza- 
tion properties of Damped Lya systems (DLAs) from the analysis of the ratio 7?.(A1 2+ /A1 + ) 
ee iV(Al 2+ )/A(Al+). We find good correlations logiV(Al + ) - log7V(Si + ) and logiV(Al + ) - 
log./V(Fe + ) that we use to indirectly estimate Af(Al + ) from iV(Si + ) and/or iV(Fe + ) measure- 
ments. In this way we determine the ratio 7?.(A1 2+ /A1 + ) for a sample of 20 DLAs. Contrary 
to common belief, the ratio can attain relatively high values (up to 0.6), suggesting that gas of 
intermediate ionization state plays an important role in DLAs. On the other hand, the lack of any 
trend between abundance ratios, such as Si/H and Si/Fe, and 7?.(A1 2+ /A1 + ) indicates that abun- 
dances are not severely influenced by ionization effects. We find a log7^(Al 2+ /Al+) - logA r (H°) 
anticorrelation that we use, in conjunction with idealized photoionization equilibrium calcula- 
tions, to constrain the ionization properties and to predict ionization corrections in DLAs. We 
consider two possible origins for the species of low and intermediate ionization state: (1) neutral 
regions devoid of Al 2+ and/or (2) partially ionized, Al 2+ -bearing regions. The log7?.(Al 2+ /Al + ) 
- log A^(H°) anticorrelation can be naturally explained in terms of a two-region model with a soft, 
stellar-type ionizing radiation field. We present abundance ionization corrections for 14 elements 
of astrophysical interest derived with different types of ionizing spectra. For most of these el- 
ements the corrections are generally below measurements errors, contrary to the predictions of 
recent models presented in the literature. We briefly discuss the potential effects of inaccuracies 
of the Al recombination rates used in the photoionization calculations. 

Subject headings: Quasars: absorption lines — Intergalactic medium — Galaxies: abundances — Radia- 
tive transfer 

1. Introduction universe. In particular, abundance studies have 
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The QSO absorption systems with hydrogen 
column densities iV(H°) > 2 x 10 20 cm~ 2 — called 
damped Lya systems (DLAs) — are believed to 
originate in intervening galaxies or protogalax- 
ies located at cosmological distances (Wolfe et 
al. 1995). High resolution spectroscopy of DLAs 
is a fundamental tool for probing the chemical 
and physical properties of the associated high-z 
galaxies and, more generally, of the high-redshift 



been performed in order to probe the nature of 
DLA galaxies and to cast light on the early stages 
of galactic evolution (Lu et al. 1996, Pettini et 
al. 1997, Prochaska & Wolfe 1999, Molaro et al. 
2000) . In abundance studies it is important to take 
into account both dust depletion and ionization ef- 
fects before deriving conclusions on the nucleosyn- 
thetic processes at work in DLA galaxies. The ef- 
fects of dust depletion have been investigated by 
several research groups (Lauroesch et al. 1996; 
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Kulkarni, Fall & Truran 1997; Vladilo 1998, here- 
after Paper I; Savaglio, Panagia & Stiavelli 2000). 
All these studies indicate that dust corrections can 
have a significant impact on our understanding of 
the chemical history of DLAs. In order to cast 
light on the intrinsic abundance trends some au- 
thors have analyzed DLAs with modest or negligi- 
ble dust content (Pettini et al. 2000, Molaro et al. 
2000) or elements which are essentially unaffected 
by dust depletions (Centurion et al. 1998). 

In the present work we investigate the ef- 
fects of ionization on the observed abundances. 
These effects are generally neglected because the 
abundances are derived from dominant ionization 
states of the elements of interest. It is well known 
from Galactic interstellar studies that the domi- 
nant ionization state in H° regions is the neutral 
one for elements with first IP > 13.6 eV (e.g. 0°, 
N°) and the singly ionized one for elements with 
first IP < 13.6 eV and second IP > 13.6 eV (e.g. 
C+ Mg+, A1+, Si+ S+ Cr+, Mn+, Fc+, and 
Zn + ). The reason for this is that the bulk of the 
H° gas is self-shielded from \\v > 13.6 eV photons, 
but transparent to hu < 13.6 eV photons. The 
first studies of ionization balance in DLAs have 
indicated that ionization corrections are negligible 
for the low ions used in abundance determina- 
tions (Viegas 1995, Lu et al. 1995, Prochaska & 
Wolfe 1996). However, the presence of Al 2+ in 
DLAs suggests that ionization corrections may be 
important. From the aluminium ionization poten- 
tials — IP(AP)=5.99 eV and IP(A1+)=18.83 eV 
- one expects Al + to be dominant and Al 2+ to 
be essentially absent in H° regions. The observa- 
tions indicate that Al 2+ is present in DLAs at the 
same radial velocity of low ions, but at different 
radial velocity of high ions such as C +3 and Si +3 
(Lu et al. 1996, Prochaska & Wolfe 1999, Wolfe & 
Prochaska 2000) . Clearly, Al 2+ is a tracer of mod- 
erately ionized gas associated to the neutral phase 
and the study of the A1 2+ /A1 + ratio must bear 
important information on the effects of ionization 
on the derived abundances. 

In a previous work we have performed detailed 
ionization calculations relevant to the ionized layer 
directly exposed to the ionizing radiation field 
(Howk & Sembach 1999; hereafter Paper II). In 
that work the photoionization calculations have 
been stopped at the point in which the local ion- 
ization fraction of neutral hydrogen climbs above 



10%. In the present work we first study the 
A1 2+ /A1 + ratio from the observational point of 
view and we then perform photoionization cal- 
culations aimed at modeling this ratio. Rather 
than analyzing individual absorbers, we study the 
general behavior of the A1 2+ /A1 + ratio in a rela- 
tively large sample of DLAs. In our computations 
we consider the possibility that the line of sight 
can intersect both neutral gas and partially ion- 
ized gas. At variance with Paper II, we also per- 
form photoionization calculations stopped at fixed 
N(tt°) values. 

2. Analysis of observational data. 

In Table 1 we present a compilation of alu- 
minium column densities in DLAs collected from 
the literature. References to the original works are 
reported in column 6. Only measurements per- 
formed on high resolution spectra obtained with 
10-m class telescopes have been included in the 
list. Most of the data have been obtained with the 
Keck HIRES spectrograph and analyzed by only 
a few authors. The collected data set is therefore 
quite homogeneous from the observational point 
of view. Nevertheless, we paid special attention 
on the reliability of the published column densi- 
ties. In particular, we revised the available data 
for rejecting cases suspected to be affected by sat- 
uration or by telluric contamination. In addition, 
we investigated the velocity distribution of the ab- 
sorption profiles. In fact the comparison of the 
column densities of A1+, Al 2+ and other low ions 
that we perform here makes sense only if these 
ions have similar velocity profiles, suggesting that 
they are physically related. 

The Al + column densities are derived from the 
transition at Ai a b 167.0 nm which is often satu- 
rated. In fact, only lower limits are reported in the 
literature for a large fraction of cases. We have, 
in addition, found three more cases of suspected 
saturation from a revision of the absorption pro- 
files. We treat also these three cases as lower lim- 
its even if they are quoted as real measurements 
in the original works (see note d in the table). 

The Al 2+ column densities are derived from 
the transitions at Ai a bl85.4nm and Ai a bl86.2nm 
which arc not affected by saturation effects. Our 
re-analysis of the velocity profiles confirms that 
Al 2+ has, in general, the same velocity distribu- 
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tion of the low ions (Wolfe & Prochaska 2000). 
However, we found some exceptions which are 
marked with the explanatory notes / and g in Ta- 
ble 1. Typically, we found evidence for an addi- 
tional Al 2+ absorption component displaced from 
the zero velocity component seen in Al + and/or 
other low ions. In most of these cases the equiv- 
alent width of this extra absorption is negligible 
compared with the total Al 2+ absorption. The 
bulk of Al 2+ originates in the same velocity range 
as Al + and the other low ions. 

For the absorber at z = 3.3901 towards QSO 
0000-263 we found evidence for telluric contami- 
nation of the Al 2+ 185.4 nm profile from the anal- 
ysis of a newly obtained spectrum of this target 
by means of UVES at VLT (Molaro et al. 2000). 
In table 1 we report the column density obtained 
from the analysis of the 186.2 nm line, which is un- 
affected by contamination. Careful inspection of 
the full set of profiles in the other DLAs of Table 1 
does not reveal evidence for telluric contamination 
since the profiles of the two lines of the doublet 
are always similar when they are both available 
(almost in all cases). 

2.1. Column-density correlation analysis 

While Al + and other singly ionized species with 
IP < 13.6 eV can be present in the inner parts of 
H° regions, the production of Al 2+ requires pho- 
tons with hv > 18.8 eV which cannot as easily 
penetrate large H° column densities. A careful 
analysis should therefore reveal a different behav- 
ior between Al 2+ and singly ionized species. Since 
the radial velocity analysis does not show evidence 
for a distinct behavior, we investigated the column 
density behavior of such species. For two species 
that are present in roughly constant proportions 
in the same region we expect to find a linear cor- 
relation with slope of unity between their loga- 
rithmic column densities. The identification of 
such a correlation would support a common ori- 
gin of the different species. The lack of such a 
correlation could be ascribed to differences in the 
abundance/depletions patterns or to differences in 
the ionization properties. By analyzing ions of the 
same elements, such as Al 2+ and Al + , there can be 
no differences in the intrinsic abundances. Study- 
ing such ions, therefore, gives information on the 
ionization state of the gas, and a slope differing 
from unity would indicate that ionization proper- 



ties do affect the observed column densities. Un- 
fortunately, only a modest number of simultaneous 
determinations of Al + and Al 2+ in DLAs arc cur- 
rently known (Table 1). To by-pass this limitation 
we compared the column densities of both Al + and 
Al 2+ ions with those of Fe + and Si + , for which a 
large number of measurements are instead avail- 
able. In practice, we performed a linear regression 
analysis of the column densities of each possible 
combination of pairs of the species Al + , Al 2+ , Si + , 
and Fe + . The results are summarized in Table 2, 
where we give for each pair the number of data 
points available, n, the correlation coefficient, r, 
the slope, to, and the intercept, q, resulting from 
the linear regression. In Fig. 1 we show the results 
for A1+ and Si+. 

One can see from Table 2 that singly ionized 
species, including Al + are well correlated with 
each other and have slope of unity within the er- 
rors 1 . This fact, together with the similarity of 
the velocity profiles, suggests that all these species 
originate in the same region. On the other hand, 
pairs of species including Al 2+ show a modest cor- 
relation coefficient and a slope significantly lower 
than 1. In addition, such pairs have larger dis- 
persions than the corresponding pairs with Al + . 
For instance er(Al 2 +,Si+) ~ 2 cr(Al+,Si+). This 
distinct behavior of Al 2+ cannot be attributed to 
differences in the abundances or dust depletions 
properties of Al, Fe and Si. If such effects were 
important, they would also tend to cancel the cor- 
relations between Al + and Fe + and between A1+ 
and Si + , which are instead clearly detected. Since 
the distinct behavior of Al 2+ is not due to intrinsic 
variations of the abundances/depletions, it must 
be due to the ionization properties of DLA clouds. 
This observational result can be interpreted in two 
ways: 

1. The Al 2+ originates in the same H° region as 

Chemical evolution and dust effects are hidden within the 
dispersion of the correlations. At first sight it is surprising 
that such effects do not introduce a large scatter. However 
it is possible to show that changes in metallicity level and 
dust content tend to compensate each other since relative 
abundances scale with the dust-to-metals ratio (Eq. 17 in 
Paper 1) and at the same time metallicity and dust con- 
tent are correlated in DLAs (Fig. 2 in Paper 1). On the 
other hand, the intrinsic spread of relative abundances due 
to chemical evolution may be modest owing to relatively 
limited range of look-back time explored in DLAs (most 
absorbers are found at z ~ 2.5). 
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the bulk of the low- ionization lines. The lo- 
cal value of the A1 2+ /A1+ ratio varies among 
DLA clouds. 

2. The Al 2+ originates in a region distinct 
from, but physically related to, the H° region 
where the bulk of the singly-ionized species 
reside. The relative column density contri- 
bution from the two regions varies among 
DLA clouds. 

Studies of the ionization properties of DLAs 
usually consider only the first of the two possi- 
bilities. Here we will also consider the second one. 
A separate Al 2+ region may in fact exist if the 
photons with hv > 18.8 eV required for producing 
Al 2+ are not available in the self-shielded parts of 
the H° region. Since the velocity structure of Al 2+ 
and of low ions are generally similar, the Al 2+ and 
H° regions must be physically connected. A cloud 
structure consisting of a partially ionized Al 2+ - 
bearing interface bordering a H° region opaque to 
hv > 13.6 eV photons satisfies the above require- 
ments. 

2.2. The A1 2 +/A1+ ratio in DLAs 

The ratio A1 2+ /A1 + can be used to constrain 
the ionization conditions and/or ionization struc- 
ture of DLAs. From the observational point of 
view we are able to measure the column density ra- 
tio 7e(Al 2 +/Al+) = 7V(Al 2 +)/7V(Al+), which gives 
information integrated along the line of sight. Di- 
rect measurements of 7?.(A1 2+ /A1+) are currently 
possible only in a few cases, mainly because reli- 
able A1+ measurements are quite rare (Table 1). 
For this reason we take advantage of the existence 
of the log7V(Al+) - log7V(Fe+) an d log7V(Al+) 
- log./V(Si + ) correlations to indirectly estimate 
A(A1 + ) from Fe + and Si + measurements. The in- 
direct determinations based on Si + measurements 
are displayed in Table 1, column 7. References to 
the original Si + measurements are the same given 
in column 6. These indirect estimates of iV(Al + ) 
yield results consistent with all the available A1+ 
data (i.e., 9 determinations and also 9 lower lim- 
its). Similar results arc obtained when Fc + is used 
to track A1+ (see more details in Section 2.4). Be- 
cause the correlation with A^(Si + ) has lower dis- 
persion, lower slope error, and lower intercept er- 
ror than the correlation with N(Fe + ) (Table 2) we 



use the Si + -bascd results given in Table 1 in the 
rest of the paper, unless otherwise specified. 

In the last column of Table 1 we list the val- 
ues of the aluminium ionization ratio in DLAs ob- 
tained by using indirect and, when possible, direct 
measurements of Al + . The ionization ratio shows 
a spread of more than one order of magnitude, 
with values ranging from 7?.(A1 2+ /A1+) ~ 0.03 up 
to ~ 0.6, and a median value ~ 0.2. Values of 
7?.(A1 2+ /A1+) as high as 0.6 are detected also in 
the sub-sample for which A^(A1 + ) is measured di- 
rectly. The presence of high fractions of doubly 
ionized aluminium suggests that A1+ is not the 
only ionization state of aluminium present in sig- 
nificant amounts in DLAs. However, it is unsafe to 
use 7?.(A1 2+ /A1+), which is a quantity integrated 
along the line of sight, as an indicator of the local 
ionization ratio. In fact, singly and doubly ionized 
aluminium may arise in distinct regions, in which 
case 7?.(A1 2+ /A1+) would reflect the relative con- 
tribution of such regions to the column densities, 
rather than the local ionization ratio. 

The 7?.(A1 2+ /A1 + ) ratio docs not show any 
trend with absorber redshift. In the interval 
2.0 < z abs 2.5, where most of the measure- 
ments are concentrated, the ratio shows the full 
spread of a factor of 20. Since metagalactic effects 
are expected to have a smooth variation with red- 
shift, the large spread in a relatively narrow red- 
shift interval suggests that the ionization fraction 
is severely influenced by local effects. The meta- 
galactic radiation field may play a role in deter- 
mining the ionization balance, but must be mod- 
ulated by some local effect. Such modulation may 
result from emission/absorption of radiation in- 
ternal to the DLAs or from large changes in cloud 
density. 

2.3. Abundances versus A1 2+ /A1+ ioniza- 
tion ratio 

The presence of high fractions of Al 2+ with 
same velocity distribution of low-ionization might 
question the reliability of abundance determi- 
nations in DLAs. If abundance measurements 
are strongly affected by ionization conditions, 
we would expect that they show some depen- 
dence on 7?.(A1 2+ /A1 + ). For instance, the models 
calculations performed in Paper II predict that 



4 



[Zn/H] and [Si/H] measurements 2 can be easily 
overestimated by one order of magnitude when 
ft(Al 2 +/Al+) > 1/10 for the relatively low ion- 
ization parameters discussed in that work. The 
collection of 7\L(A1 2+ /A1 + ) values presented in 
Table 1 allows us to probe if such effects are 
indeed present, at least for elements commonly 
measured in DLAs. We performed a linear re- 
gression analysis of [Si/H] and [Si/Fe] against 
logft(Al 2 +/Al+). The [Si/H] analysis yields cor- 
relation coefficient r = 0.29, dispersion a — 0.46, 
and slope m = 0.37 ± 0.30 (18 data points). The 
predicted increase of [Si/H] at high 7£(A1 2+ /A1 + ), 
if present, is not significant from the statistical 
point of view. Similar results are obtained from 
the study of [Si/Fe], a ratio which is expected 
to be moderately sensitive to ionization effects 
(Paper II). The [Si/Fe] analysis yields r = 0.05, 
cr = 0.20, and m = -0.03 ± 0.13 (19 data points). 
In this case the null result is even more clear. The 
lower dispersion than in the case of the Si/H ra- 
tio is probably due to the fact that variations of 
the metallicity and of the dust-to-gas ratio tend 
to cancel when we consider a relative abundance 
such as Si/Fe. These results suggest that ioniza- 
tion corrections arc not severe, or at least that 
they do not show a strong, obvious dependence on 
ft(Al 2 +/Al+). 

2.4. A1 2+ /A1+ ratio versus H° column den- 
sity 

In Fig. 2 we plot ft(Al 2 +/Al+) versus the ab- 
sorber H° column density. The data points show a 
general decrease of log 7?.(A1 2+ /A1 + ) with increas- 
ing logA(H°). A linear regression analysis yields 
an anti-correlation with slope m = —0.81 ± 0.15 
and intercept q = 16.0 ±3.1. Even if the cor- 
relation coefficient is not very high (Pearson's 
r = —0.80), the probability of the null hypoth- 
esis (of no correlation) is 5.9 x 10 -5 . The highest 
value of the ratio, ^(A1 2 +/A1+) ~ 0.6, is found 
at A(H°) ~ 10 20 - 3 cm" 2 , while the lowest value, 
ft(Al 2 +/Al+) ~ 0.03, is found at A(H°) ~ 10 21 ' 4 
cm . 

The above result is based in large part on 
indirect Al + determinations obtained from Si + 
measurements. However, the log 7?.(A1 2+ /A1 + ) 

2 We adopt the usual convention [X/Y] = log[ j^jyj ] — 
log(X/Y) Q . 



logA(H°) anti-correlation is confirmed when Fe + 
is used to infer Al + column densities. In fact, 
the slope and intercept derived in this case - 
m = —0.78 ± 0.16 and q = 15.2 ± 3.3 — are very 
similar to the above given values, and the proba- 
bility of the null hypothesis equals 8.6 x 10~ 5 . 

The lack of DLAs at the bottom-left and top- 
right corners of Fig. 2 could be due, in princi- 
ple, to the impossibility of detecting Al + and/or 
Al 2+ below their observational limits. However, 
we have verified that this is not the case. In fact, 
the strength of the Al + 167.0 nm transition and of 
the numerous Si + and Fe + transitions used for the 
indirect Al + determinations guarantee that Al + 
would be easily detected down to much lower col- 
umn densities than observed. In the case of Al 2+ , 
the detection limit with Keck+HIRES is ~ 10 12 
cm~ 2 . In the sample of Table 1, only one case out 
of 21 measurements has such a low value. Of the 
remaining A(A1 2+ ) data, 9 are 0.5 dex above the 
detection limit and 11 are 1 dex above the limit. 
We conclude that the anti-correlation is not an 
artifact induced by observational limitations. 

In addition, we do not have reasons to believe 
that the observed trend results from a selection 
bias. The most common bias considered to affect 
the population of DLAs is dust obscuration of the 
background QSO (Pei et al. 1991). In fact, this 
effect may be responsible for the apparent anti- 
correlation between [Zn/H] and iV(H ) found by 
Boisse et al. (1998), if lines of sight of high metal- 
licity (dust content) and high column density in- 
deed obscure the background QSO. However, the 
same effect should not be relevant for the anti- 
correlation reported here given the fact that the 
ratio 7£(A1 2+ /A1 + ) is independent of metallicity. 

The existence of an empirical anti-correlation 
log N(B°) - logft(Al 2+ /Al+) implies that the neu- 
tral hydrogen column density can be used as an 
indirect estimator of the ionization state of the 
gas. In this way one can search for ionization ef- 
fects in DLAs without available Al measurements. 
Given the importance of Zn in DLA studies we 
searched for a possible dependence of [Cr/Zn] and 
[Fe/Zn] on log A(H°). A linear regression through 
the 26 [Cr/Zn] measurements available yields a 
slope m — —0.03 ± 0.11 and a correlation coef- 
ficient r = 0.06. Very similar results are obtained 
from the analysis of the 22 [Fe/Zn] available de- 
terminations: m = -0.04 ± 0.13 and r = 0.07. In 
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addition, there is no trace of an increased disper- 
sion at low iV(H°) for any of the two ratios. These 
empirical results suggest that Zn ionization effects 
are probably negligible in DLAs. 

If A1 2+ /A1 + is relatively constant throughout 
the cloud, the 7?.(A1 2+ /A1+) ratio should be ap- 
proximately constant with AT(H°). The decrease 
of -ft(Al 2 +/Al+) with N(U°) can be interpreted in 
two ways, depending on which of the two possible 
origins of Al 2+ mentioned at the end of Section 
2.1 is more appropriate. (1) If Al 2+ is co-spatial 
with the rest of low- ionization species, the anticor- 
relation is consistent with a reduction of the gas 
ionization level with increasing H° self-shielding. 
(2) If Al 2+ originates in a partially ionized inter- 
face bordering a neutral region, then we expect 
the column density of the ionized interface and of 
the neutral region to be unrelated and therefore 
A(A1 2 +) to be independent of JV(H°). On the 
other hand, we expect A1+ and other species of 
low ionization to scale with H° owing to a com- 
mon origin in the neutral region. The decrease of 
ft(Al 2 +/Al+) with N(W) is naturally explained 
from the fact that A^(A1 2+ ) does not scale with 
iV(H°) while N(A1+) does. 

3. Model calculations 

On the basis of the observational results dis- 
cussed above, we consider two possible origins 
for species of low and intermediate ionization: 
(1) regions completely opaque to photons with 
hv > 13.6 eV and/or (2) regions partially trans- 
parent to ionizing photons. Neutral and singly 
ionized species can arise in both regions, but Al 2+ 
can only be present in typc-2 regions. The gen- 
eral idea is that type-2 regions are the photoion- 
ized envelopes of type-1 regions. The similarity 
of the velocity profiles of Al 2+ and singly ionized 
species is naturally explained in this way. On the 
other hand, radial- velocity studies show a general 
misalignment of C +3 and Si +3 profiles relative to 
Al 2+ profiles, suggesting that the bulk of high- 
ionization species originates elsewhere (Wolfe & 
Prochaska 2000). We expect therefore that type-2 
regions are mildly ionized and have low Si +3 /Si + 
fraction. 

In the Appendix we derive the expressions that 
allow us to compute ionization ratios and abun- 
dance ionization corrections in the framework of 



the two-region model. These expressions do not 
depend on the metallicity, provided type-1 and 
type-2 regions have equal metallicity in any given 
DLA system. The relative contribution of the two 
regions along the line of sight is specified by the pa- 
rameter iVi /N 2 , which is the ratio of the total col- 
umn densities in region 1 and 2, respectively. The 
possibility that Al 2+ and the species of low ioniza- 
tion originate together in a single layer is included 
in our treatment as the special case Ni/N% = 0. In 
fact, the single layer should be a partially ionized 
region of type 2 given the presence of significant 
fractions of Al 2+ . 

In order to model the ionization properties of 
the gas we assume that the DLA regions under 
study are embedded in a ionizing radiation field 
of given spectrum and intensity. In type-1 regions 
we only consider the contribution of species which 
are dominant ionization states in HI regions. In 
type-2 regions we compute the ionization fractions 
by means of photoionization equilibrium calcula- 
tions. For this purpose we used the CLOUDY 
code (v90.04; Fcrland 1996, Ferland et al. 1998) 
assuming plane-parallel geometry with ionizing ra- 
diation incident on one side. We consider two 
possible types of radiation fields: a hard, QSO- 
dominated spectrum representative of the radia- 
tion field external to the DLAs at z — 2 (Haardt & 
Madau 1996; Madau, Haardt & Rees 1999) and a 
soft, stellar-type spectrum (T e ff =33,000 K; Kurucz 
1991) representative of the internal radiation field 
or of an external field dominated by starlight from 
galaxies (e.g. Steidel, Pettini & Adelberger 2000). 
In both cases the intensity of the field is speci- 
fied by the ionization parameter U = $(H)/cnn , 
where ^(H) is the total surface flux of ionizing 
photons (cm~ 2 s _1 ) and njj the hydrogen par- 
ticle density (cm -3 ). We refer to Paper II for 
more details on the model assumptions and on the 
adopted radiation fields. In the rest of this section 
we present model predictions of the ionization ra- 
tios 7?.(A1 2+ /A1+) and of the abundance ionization 
corrections. 

3.1. Aluminium ionization ratio 

3.1.1. Soft-continuum, two-regions model (S2 
model) 

The photoionization calculations indicate that 
with the soft, stellar-type spectrum the H° column 
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density of the ionized layer, A^H ), cannot attain 
the high values typical of DLAs for a wide range 
of U. In fact, we find iV 2 (H°) < 3.2 x 10 18 cm- 2 
for —5.4 < logf/ < —1.0. Therefore, we assume 
that a neutral region of type 1 is present along 
the line of sight, in addition to the Al 2+ -bearing, 
type-2 region. We refer to this soft-continuum, 
two-regions model as S2. 

Because the ionized region does not yield a sig- 
nificant contribution to the neutral hydrogen col- 
umn density and because x(H°) = 1 in the neutral 
region, we take Aq(H) = Aq(H°)= N(tt°) as the 
total column density of the type-1 region, where 
N(H°) covers the interval 10 20 2 cuT 2 < N(R°) < 
10 22 cm -2 , representative of DLAs. For each given 
A^i (H) we let the ionization parameter U vary and 
we use CLOUDY to determine the column den- 
sity of the type-2 region, N 2 (H.), and the mean 
ionization fractions of the species of interest. We 
then estimate ft(Al 2+ /Al+) by means of Eq. 7, 
where we take N 1 /N 2 = Ni (H) / N 2 (H) and we 
insert the Al + and Al 2+ ionization fractions. In 
this way we derive model curves of 7^(A1 2+ /A1+) 
versus N(H°) that can be compared with the ob- 
served data points. 

In Fig. 3 we show the results of these calcu- 
lations for constant values of the ionization pa- 
rameter. The solutions calculated at U = 10~ 2,2 
(thick curve) are consistent, within the statisti- 
cal errors, with the linear regression of the ob- 
servational data points (dashed-dotted line). The 
corresponding total hydrogen column density of 
the type-2 layer predicted by the model is A^(H) 
= 8 x 10 20 cm -2 . The solutions calculated at 
t^min = 1CP 2 ' 6 and /7 max = 1CP 1 ' 7 (lower and up- 
per thin curves, respectively) bracket the obser- 
vational data points. For these lower and upper 
envelope curves we find A^(H) = 3.3 x 10 20 and 
2.3 x 10 21 cm -2 , respectively. The total column 
densities of type-1 and type-2 regions are roughly 
comparable (N1/N2 ~ 1). It is remarkable that 
the logAT(H°) - log^(Al 2 +/Al+) anticorrelation 
can be easily modeled with very simple assump- 
tions and with a single interval of U values in the 
framework of the S2 model. 

3.1.2. Hard- continuum, single-region model (HI 
model) 

With the QSO-dominated field the H° column 
density of the ionized layer can attain the high val- 



ues typical of DLAs if U is sufficiently high. As 
a consequence, the presence of a of type-1 region 
is not required in order to explain the observed 
AT(H°) values. In this case we assume that all 
the low and intermediate ionization species orig- 
inate in a single, partially ionized layer. We re- 
fer to this hard-continuum, one-region model as 
HI. We stopped the calculations at a series of 
A^H ) values representative of DLAs, rather than 
at a given threshold of ionization fraction. Our 
treatment of this case is therefore quite similar 
to those of Lu ct al. (1995) and Prochaska & 
Wolfe (1996). At variance with the S2 model, in 
the HI model there is residual ionized gas, and 
hence Al 2+ , throughout the interior of the cloud 
due to the strong high-energy tail of the extra- 
galactic Haardt & Madau spectrum. In Fig. 4 
we show the resulting 7£(A1 2+ /A1 + ) ratios plot- 
ted versus AT(H°) at constant values of U. We 
estimated 7\L(A1 2+ /A1+) from Eq. 7 by taking 
Ni/N 2 = 0. When U is too low there are diffi- 
culties with the thermal solution in the CLOUDY 
computation since the temperature falls well below 
1000 K before the total N(H°) is reached. This ex- 
plains the lack of solutions below log U = —4.8 in 
Fig. 4. On the other hand, some of the solutions 
tend to overproduce Si +3 when U is too high. As 
an example, in Fig. 4 we indicate with empty cir- 
cles the solutions for which ft(Si 3+ /Si+) > -0.5 
dex. The intrinsic 7?.(Si 3+ /Si + ) is probably much 
lower than this conservative limit given the very 
different velocity structure of Si 3+ and Si + . 

With the HI model the curves predicted at 
constant U have very different - in some cases 
opposite - slope from that of the anticorrelation 
logAT(H°) - log"K(Al 2 +/Al+) (Fig. 4). When all 
the U values are considered, the solutions tend to 
fill the plane log A^(H°) - logft(Al 2 +/Al+) without 
any preference for the regions populated by obser- 
vational data. To reproduce the anti-correlation it 
is necessary to impose very specific constraints to 
the input parameters. For a given A^(H°), there 
is an allowed interval of U values such that the 
predicted log7?.(Al 2+ /Al + ) ratios overlap the ob- 
served ones. For instance, at log N(H°) = 20.4 we 
must require that —4.8 < log U < —3.6 in order to 
obtain ratios that lie within ±la of the regression 
to the observed data points (Fig. 4). In addi- 
tion, the ionization parameter must decrease, on 
the average, while N(H°) increases. For instance, 
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we must require that U oc 7V(H°) 15 in order to 
reproduce the slope of the observed anticorrela- 
tion in the range 20.4 < logN(H°) < 20.8. This 
is at variance with the results of the S2 model, 
for which all the observed data points are easily 
matched by adopting a single interval of U values. 

The use of one-side illuminated clouds is in- 
herent to the design of CLOUDY, but may be 
slightly inappropriate for the external, hard ion- 
izing spectrum case. According to Prochaska & 
Wolfe (1996) 1-side calculations tend to give a 
lower degree of ionization than 2-side calculations. 
We have estimated this effect by doubling the col- 
umn densities of the 1-side calculations. With this 
type of estimate we do not find significant changes 
in the results of the present work, including the 
abundance corrections discussed below. 

3.2. Abundance ionization corrections 

Abundance determinations in DLAs are based 
on column-density measurements of species which 
are dominant ionization states in H° regions. In 
the Appendix we define the ionization correction 
terms that allow us to recover the intrinsic abun- 
dances starting from this type of measurements. 
In order to estimate these correction terms we used 
the same sets of input parameters — ionizing spec- 
trum, U and Ni/N 2 — that allow us to repro- 
duce the observed {A(H°),^(A1 2+ /A1+)} distri- 
bution. In practice, with such parameters we de- 
termine the mean ionization fractions ^(X 1 ) with 
CLOUDY and the correction terms from Eqs. 5 
and 6. Thanks to this procedure the input param- 
eters are constrained by the requirement to model 
the A1 2+ /A1 + observations. An additional advan- 
tage is that the ionization effects are estimated as 
a function of N(H°), which is a measurable quan- 
tity. In this way, the range of predicted ionization 
corrections is reduced once we know AT(H°) for in- 
dividual DLAs. The resulting corrections for abso- 
lute abundances are shown in Figs. 5 through 16. 
In each figure we display the corrections predicted 
by the S2 and HI models discussed above. Cor- 
rections for N/H and O/H are not shown in the 
figures because they are always below abundance 
measurement errors ( < 0.05 dex). 

The corrections estimated with the S2 model 
are shown in Figs. 5 through 16 (absolute abun- 
dances) and Table 3 (relative abundances). In the 
S2 model the species used for abundance measure- 



ments arise mostly in type-1 neutral regions, with 
an additional contribution from type-2 ionized re- 
gions. It is this latter contribution that affects 
the measured abundances. The correction terms 
predicted by this model are generally small owing 
to the high U values that we require to match the 
7e(Al 2 +/Al+) versus A(H°) data. In fact, when U 
is high, the species used for abundance measure- 
ments tend to vanish in the ionized region since 
they move to a higher state of ionization. Ioniza- 
tion corrections tend to decrease in absolute value 
with increasing N(H°). This is due to the fact 
that the relative contribution of the ionized layer 
becomes less important at high N(H°). 

The predictions for the HI model have been 
derived by considering only the solutions which 
match the empirical log Y(H°) - log ft(Al 2 +/Al + ) 
anticorrelation. The empty circles in Figs. 5 
through 16 have been calculated using the solu- 
tions found at the intersection between the linear 
regression to the data points and the curves at 
constant ionization parameter log{7 = —4.2 and 
—4.8 shown in Fig. 4. Owing to the difficulty 
of finding solutions at lower U and high Y(H°) wc 
have not calculated corrections at log N(H°) > 21. 
It is clear, however, that the correction terms be- 
come negligible at high log A(H°). The reason for 
this is that the corrections decrease (in absolute 
value) with decreasing U and, at the same time, 
U oc N(H°)~ 1,5 . The predicted corrections are 
generally small because the match with the ob- 
served 7?.(Al +2 /Al + ) ratios is found at low values 
off/. 

The correction terms can be negative or posi- 
tive depending on the model adopted and on the 
species considered. The corrections for absolute 
abundances X/H are negative in the S2 model. 
In this model the dominant species used for the 
measurements, X ld and H°, mostly arises in the 
type-1 region unaffected by ionization; the type-2 
layer gives an extra contribution which enhances 
the X ld /H° ratio and a negative correction is re- 
quired to recover the intrinsic abundance. In the 
HI model the corrections for absolute abundances 
can be negative or positive since all the species 
arise in a single layer in which the X ld /H° ratio 
can be larger or smaller than the intrinsic abun- 
dance. 
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3.3. Accuracy of Al atomic data 

The above calculations rely on the accuracy of 
the input atomic data. Here we are particularly in- 
terested in the accuracy of Al data since our mod- 
els are constrained by the capability of matching 
the observed log A(H°)-log:K(Al 2 +/Al+) anticor- 
relation. The photoionization cross-sections cur- 
rently used in CLOUDY calculations are generally 
accurate within ~ 10%, the A1+ cross section hav- 
ing a regular level of accuracy (Verner et al. 1996; 
Ferlandetal. 1998). Radiative recombination rate 
coefficients can be obtained with an accuracy bet- 
ter than ~ 15% (Ferland et al. 1998). However, 
the recombination process can be dominated by 
dielectronic recombination (DR), which is a far 
more uncertain mechanism. Low-temperature DR 
rates, which are critical for determining the ion- 
ization balance in photoionization equilibrium, are 
lacking for many elements (Ferland et al. 1998). 
Luckily, such coefficients have been calculated for 
Al. However, given the theoretical and experi- 
mental uncertainties it is possible that the A1+ 
DR rate may be overestimated (Nussbaumer & 
Storey 1986). If the Al + recombination rate is 
overestimated, the predicted 7?.(A1 2+ /A1+) ratio 
is underestimated. An effect of this type has been 
reported in a photoionization study of a Lyman 
Limit system, in which the models that give a 
good fit to other species are not able to repro- 
duce the relatively high 7?.(A1 2+ /A1 + ) ratio ob- 
served (D'Odorico & Petitjean 2001). In order to 
test the consequences of an effect of this type we 
artificially increased the ratio 1Z(A\ 2+ / A\ + ) calcu- 
lated at any given value of U. The results that we 
found can be summarized as follows. In the S2 
model (i) the observed anticorrelation is matched 
at lower U values; (ii) the column density N 2 of the 
ionized layer, which scales with U, becomes lower; 
(iii) the ionization corrections, which scale with 
N 2 , become lower; (iv) the Al correction becomes 
even lower (in absolute value) owing to the reduced 
contribution of Al + from the ionized layer. In the 
HI model (i) the solutions that match the observed 
anticorrelation are shifted to lower U values; (ii) 
the ionization corrections, which in this case scale 
with U, become lower. In summary, if the Al re- 
combination rate is too high, then the ionization 
parameter U and the abundance corrections calcu- 
lated above should be reduced; in particular, the 
Al abundance corrections of model S2 would be 



more in line with those of the other elements. 

3.4. Implications for the indirect esti- 
mates of N(Al+) 

The large values of Al correction terms that we 
find are somewhat surprising given the fact that 
most A1+ column densities have been indirectly es- 
timated from Si + column densities. 3 In fact, the 
log iV(Al + )-log ./V(Si + ) correlation could, in prin- 
ciple, be destroyed by the large ionization effects 
predicted. Since this is not the case, we must un- 
derstand why. One possible reason is that, ow- 
ing to the uncertainty of Al + recombination coef- 
ficients, the Al correction terms may be overesti- 
mated, as we have discussed above. Another pos- 
sibility is that the Al correction terms are correct, 
but the correlation is not destroyed because we are 
considering a limited interval of N(H°). In fact, 
in the column-density range of the DLAs used for 
deriving the logiV(Al + )-logiV(Si + ) correlation — 
i.e., 20.3 < log N(H°) < 20.7 — the Al/Si correc- 
tions terms are large, but nevertheless they only 
show a modest variation (f» ±0.1 dex). As a con- 
sequence, the ionization effects may significantly 
change the intercept of the logiV(Al + )-log A^(Si + ) 
correlation, but should not change significantly the 
dispersion and the slope. 

The A(A1 + ) values inferred from A(Si + ) may 
contain a systematic error when logiV(H°) > 20.7 
owing to the variation of C (Al/Si) with A(H°). 
If the Al correction terms are overestimated, 
the effect is probably small. Otherwise, we can 
quantify this error from the predicted variation 
of logC(Al/Si) between logA(H°) ~ 20.5 and 
logA(H°) ~ 21.5. We find that N(Al+) may 
be overestimated by ~ 0.2/0.3 dex in model S2 
or underestimated by a similar amount in model 
HI. These values are generally within the error 
bars quoted in the last column of Table 1. In 
any case, these errors would not affect the results 
of the present work. In the case of model S2, in 
order to compensate for such errors we should de- 
crease A(A1+) and hence increase 7?.(A1 2+ /A1+) 
by ~ 0.2/0.3 dex at logA(H°) > 20.7. One can 
see from Fig. 3 that such increase would even 
improve the agreement between the S2 model and 
the observations. In the case of the HI model, 

3 This discussion could be equally applied to the indirect 
estimates based on Fc+ column densities. 
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we should instead decrease 7\L(A1 2+ /A1 + ). This 
would imply that the interval (£7 m in,£7 ma x) of al- 
lowed solutions should be slightly shifted to lower 
U values at high logiV(H°). 

4. Implications for abundance studies 

We briefly discuss how the ionization correction 
terms presented here can affect studies of DLA 
abundances. Unless differently specified, the re- 
sults summarized below are relative to both types 
of models considered in this work. 

The correction terms for N and O are generally 
negligible. This is not surprising since the ioniza- 
tion fractions of H°, N° and 0° are held together 
because of the strong charge exchange reactions 
(e.g. Sembach et al. 2000). We note, however, 
that cosmic rays ionization can affect the N°/0° 
ratio (Viegas 1995). The [N/Si] ratio is slightly 
underestimated (Table 3). The correction term 
for [N/S] can be positive or negative, depending 
on the adopted model. In any case the effect is 
generally small, with |logC(N/S)| < +0.2 dex 
for N(U°) > 20.3. These results indicate that the 
large scatter of [N/Si] and [N/S] abundances found 
by Lu, Sargent & (1998) and by Centurion et al. 
(1998) is a genuine nucleosynthetic effect. Claims 
that nitrogen abundances are severely affected by 
ionization effects (Izotov & Thuan 1999) are not 
supported by our study. 

The correction terms for Al may be quite large, 
but with different signs depending on the model 
adopted (Fig. 6). Measurements of the Al/Fe and 
Al/Si ratios in DLAs yield [Al/Fe] sa 0, with a 
scatter of ~ 0.3 dex, and [Al/Si] w —0.4 dex, with 
values between —0.1 and —0.6 dex (Prochaska & 
Wolfe 1999). Local ISM studies indicate that Al 
and Fe have similar depletion (Barker et al. 1984; 
Howk, Savage & Fabian 1999). The [Al/Fe] ra- 
tios should not be affected by depletion, while the 
[Al/Si] ones should be affected as much as the 
[Fe/Si] ratios, which are underestimated by ~ 0.3 
dex in DLAs (Paper I). So we expect [Al/Fe] w 
and [Al/Si] w —0.1 dex after correcting for dust. 
For the DLAs with available Al + measurements 
we predict log C (Al/Fe) « logC(Al/Si) w -0.4 dex 
(S2 model) or « +0.2 dex (HI model). Therefore, 
after dust and ionization corrections are applied, 
the [Al/Fe] and [Al/Si] ratios are below solar in 
the S2 model, but somewhat enhanced in the HI 



model. As we mentioned above, these ionization 
effects may be weaker if the Al + recombination 
rate is overestimated. 

Silicon and sulphur are used as a tracers of al- 
pha elements given the difficulty of measuring oxy- 
gen in DLAs. The [Si/H] correction term is ~ —0.1 
dex at log A(H°) ~ 20.3 (Fig. 7) and silicon abun- 
dances can be accordingly overestimated. How- 
ever, the effect is not strong and this may explain 
why [Si/H] measurements in DLAs do not show 
a trend with 7\L(A1 2+ /A1 + ), as discussed in Sec- 
tion 2.3. The [S/H] corrections terms are negative 
for the S2 model and positive for the HI model, in 
both cases being < 0.2 dex in absolute value (Fig. 
9). 

The [Si/Fe] ratio, which is used as a proxy of the 
a/Fe ratio, shows an enhancement of ~ +0.3/+0.5 
dex in DLAs which has been interpreted as an in- 
trinsic nucleosynthetic effect (e.g. Lu et al. 1996) 
or as a differential dust depletion (e.g. Paper I). 
The present study indicates that [Si/Fe] is almost 
unaffected by ionization corrections, being possi- 
bly overestimated by ~ 0.1 dex at low A(H°) in 
the S2 model. Owing to the negligible dust deple- 
tion of both sulphur and zinc, the [S/Zn] ratio has 
been used as a dust-free [a/Fe] indicator in DLAs 
(Centurion et al. 2000; see however Prochaska et 
al. 2000). The correction terms for the [S/Zn] ra- 
tio are « 0.2 dex in absolute value for the DLAs 
with available S measurements. Considering that 
such corrections are dominated by the zinc contri- 
bution which may be overestimated (see below), 
we conclude that the [S/Zn] results presented by 
Centurion et al.( 2000) are modestly affected by 
ionization effects. Local interstellar studies yield 
a S/Zn ratio approximately solar (Howk, Savage 
& Sembach 1999) suggesting that ionization cor- 
rections are unimportant. 

Ionization effects tend to lower the measured 
Ar abundances in the local interstellar medium, 
at least in lines of sight with log7V(H°) < 20.0 
(Sofia & Jenkins 1998, Jenkins et al. 2000). In 
DLAs the Ar corrections are negligible in the case 
of the S2 model; however, Ar may be severely un- 
derestimated in the HI model when the H° column 
density is low (Fig. 10). 

Corrections for Cr and Fe are < 0.05 dex in 
absolute value (Figs. 12 and 14) and those for Ti, 
Mn, and Ni arc < 0.1 dex (Fig. 11, 13, and 15). 
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Therefore, deviations from solar ratios observed 
for pairs of iron-peak elements such as [Cr/Fe] and 
[Mn/Fe] (Lu et al. 1996; Prochaska & Wolfe 1999) 
must be due cither to dust depletion or to nucle- 
osynthetic effects since ionization effects are ex- 
cluded. 

The predicted corrections for Zn are relatively 
large, with opposite signs depending on the model 
adopted (Fig. 16). Owing to the predicted vari- 
ation of the Zn correction terms with log7V(H°) 
we would expect to find some trend between the 
[Cr/Zn] and [Fe/Zn] ratios and logA(H°). The 
lack of any trend (see Section 2.4) suggests that 
the zinc correction terms may be overestimated. 
This could be the case since the Zn recombination 
coefficients and ionization cross sections are rather 
uncertain (Howk, Savage & Fabian 1999). 

5. Summary and conclusions 

We have analyzed column density measure- 
ments of Al 2+ and of singly-ionized species avail- 
able in literature to cast light on the proper- 
ties of low-ionization regions in DLA systems. 
We have found that logiV(Al + ) is well corre- 
lated with log A(Si + ) and log N(Fe + ) and we have 
used this result to estimate A(A1 + ) and the ratio 
ft(Al 2 +/Al+) = 7V(Al 2 +)/A(Al+) for a sample of 
20 DLAs. The ratio can attain relatively high val- 
ues, up to ~ 0.6, with a median value of 0.2. This 
result is contrary to the common belief that the 
fraction of Al 2+ is generally small in DLAs. In 
the redshift interval 2.0 < z a b s < 2.5, where most 
of the measurements are concentrated, the ratio 
shows the full spread of a factor of 20. There- 
fore, the ratio must be influenced by local effects 
not dependent on z a b s . Local absorption and/or 
local radiation fields probably play an important 
role in determining the ionization properties inside 
DLAs. 

The presence of high fractions of Al 2+ with 
same velocity distribution of low-ionization may 
question the reliability of abundance determina- 
tions in DLAs. We have investigated the behavior 
of [Si/H] and [Si/Fe] abundance ratios in order to 
put in evidence possible effects of ionization. How- 
ever, we do not find any trend between these abun- 
dance ratios and 7^(A1 2+ /A1 + ). We have consid- 
ered the possibility that Al 2+ originates in a region 
different from the one where species of low ioniza- 



tion arise. In this case, the ratio 7£(A1 2+ /A1 + ) 
would measure the relative contribution of differ- 
ent regions, rather than the intrinsic degree of ion- 
ization of a single region. From a regression anal- 
ysis of the logarithmic column densities we find 
some evidence for a distinct behavior of Al 2+ . In 
fact, while any pair of species including Al + , Si + 
and Fe + yields correlations with slope unity and 
low dispersion, pairs including Al 2+ yield corre- 
lations with lower slopes larger dispersions. One 
possible explanation of this observational result is 
that, indeed, Al 2+ and singly-ionized species orig- 
inate in two different regions. However, the two 
regions must be physically connected since the ve- 
locity profiles of Al 2+ and of low ions are very 
similar. 

We have also identified the existence of an 
anticorrelation between log7\L(Al 2+ /Al + ) and 
logA(H°). The anticorrelation appears to be an 
intrinsic property of DLAs not induced by obser- 
vational bias, at least as far as the detection limits 
of Al + and Al 2+ are concerned. We have used such 
anticorrelation, in conjunction with photoioniza- 
tion equilibrium calculations, to constrain the ion- 
ization parameter in DLAs and hence the abun- 
dance ionization corrections. We have proposed 
that low-ionization species in DLAs may arise in 
two types of regions: (1) an H° region opaque to 
hv > 13.6 photons and/or (2) a partially ionized, 
Al 2+ -bearing interface with small/negligible frac- 
tions of high ions such as Si +3 . We have considered 
two types of ionizing continuum: a soft, stellar- 
type (T e g=33,000 K) and a hard, QSO-dominated 
type at z w 2. 

We have succesfully reproduced the observed 
logft(Al 2+ /Al+) -log A(H°) anticorrelation by 
means of a soft- continuum, two-region model (S2 
model) with ionization parameter in the range 
10~ 2 - 6 < U < lO" 17 . In this model most of 
the neutral hydrogen and low ionization species 
originates in the neutral region of type 1. How- 
ever, the total hydrogen column density of the 
partially ionized, type-2 region is relatively high, 
3 x 10 20 cm- 2 < A 2 (H) < 2 x 10 21 c m - 2 . At a 
given value of U, A 2 (H) is fixed by the photoion- 
ization calculation and the contribution of the 
type-2 region decreases as A(H°) increases. Be- 
cause Al 2+ originates only in the type-2 region, 
the anticorrelation is naturally explained. 

We also tried to reproduce the observed log7?.(Al 2+ /Al + ) 
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-log A(H°) anticorrelation by means of a hard- 
continuum, one- region model (HI model). In 
this model all the species of low and interme- 
diate ionization state are assumed to arise in a 
single layer. In order to reproduce the decrease of 
log7?.(Al 2+ /Al + ) with the correct slope it is nec- 
essary to assume that U must decrease with a law 
of the type U oc Ar(H°) -1 ' 5 . The anticorrelation 
might be due to the decrease of the mean ioniza- 
tion level with increasing self-shielding by neutral 
hydrogen. The typical values of the ionization 
parameter at logN(H°) ~ 20.8 is \ogU ~ -4.8. 

We have estimated abundance ionization cor- 
rections for 14 elements of astrophysical interest 
both with the S2 and the HI model. In both 
cases we used the same sets of parameters that 
allow us to reproduce the anticorrelation between 
logft(Al 2+ /Al+) and log N(B°). Ionization cor- 
rections can be negative or positive depending on 
the model adopted and on the species considered. 
In any case ionization corrections tend to become 
smaller in absolute value as iV(H°) increases. 

Ionization corrections are small in both mod- 
els, but for different reasons. In the S2 model 
corrections are small because the species used for 
abundance measurements tend to shift to a higher 
ionization state in the type-2 region owing to the 
high value of U that we find. In the HI model cor- 
rections are small because the single region where 
low-ionization species, together with Al 2+ , are lo- 
cated has a low level of ionization. 

The correction terms for the absolute abun- 
dances of N, O, Ti, Cr, Mn, Fe, and Ni are gener- 
ally below measurement errors (w 0.05/0.1 dex), 
independently of the adopted ionizing spectrum. 
Therefore the deviations from solar ratio observed 
in some pairs of iron-peak elements, such as the 
[Mn/Fe] or [Cr/Fe] ratios (Lu et al. 1996, Pet- 
tini et al. 2000) are not induced from ioniza- 
tion effects. The Ar/H correction term is neg- 
ligible for the S2 model but may be significant 
for the HI model, in which case the measured 
Ar°/H° value would underestimate the intrinsic 
Ar/H abundance. 

The ionization corrections for Mg, Si, P, and S 
can attain values somewhat higher than the mea- 
surement errors. The [N/Si] and [N/S] ratios are 
modestly affected by ionization effects. As a con- 
sequence, the considerable [N/S] and [N/Si] scat- 
ter observed at a given metallicity (Centurion et 



al. 1998; Lu, Sargent & Barlow 1998) is a gen- 
uine characteristic of DLAs. The [Si/Fe] ratio, a 
typical indicator of the a/Fe-peak ratio, may be 
overestimated by w 0.1 dex at low iV(H°) if ion- 
ization corrections are not applied. The [S/Zn] 
ratio might be more sensitive to ionization effects, 
but the result is uncertain since Zn corrections are 
probably inaccurate. 

The Al corrections can be relatively large. They 
can be negative or positive, depending whether 
we adopt the S2 or the HI model, respectively. 
The [Al/Fc] and [Al/Si] ratios corrected for dust 
and ionization effects are below the solar value in 
the S2 model and somewhat enhanced in the HI 
model. These effects, however, are less marked if 
the Al recombination rate is overestimated. 

The Zn corrections are apparently large, but 
these results may be inaccurate owing to the 
uncertainties of Zn atomic parameters. From 
an analysis of [Cr/Zn] and [Fe/Zn] data versus 
log7V(H°) we have provided evidence that zinc 
ionization corrections are likely to be overesti- 
mated. Therefore studies of DLAs metallicity 
based on [Zn/H] data (Pettini et al. 1997; Vladilo 
et al. 2000) are unlikely to be significantly affected 
by ionization effects. 

We have investigated the stability of the above 
results in light of possible inaccuracies of Al atomic 
parameters. In particular we have considered the 
possibility, consistent with available data, that the 
Al + dielectronic recombination rate may be over- 
estimated. We find that in this case the ionization 
parameter U and the abundance corrections would 
be lower both in the S2 and in the HI model. 

The ionization corrections presented here are 
significantly smaller than the ones predicted by 
Izotov, Schaerer & Charbonnel (2000), who have 
also considered a two-region model of DLA gas. 
However, these authors assume that the neutral 
region has much lower metallicity than the ion- 
ized one, a strong assumption for which there is 
little observational support (see Levshakov, Kegel 
& Agafonova, 2000). With such assumption the 
metal absorptions originate essentially in the ion- 
ized region and the predicted ionization effects are 
obviously more enhanced than in our case. 

Future studies of ionization properties in DLAs 
would benefit from measurements of ionic ra- 
tios other than A1 2+ /A1 + . One possibility is the 
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Fe 2+ /Fe + given the presence of the 112.2 nm tran- 
sition of Fe 2+ which could be observed in selected 
cases. From the point of view of the atomic data 
it is important to better understand the actual 
accuracy of atomic parameters, and in particular 
those of Al and Zn. 
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the Johns Hopkins University. We thank the ref- 
eree for suggestions that have significantly im- 
proved the quality of this work. 

6. Appendix: ionization corrections and 
ionization ratios in the two-region model 

We assume that low ionization species in DLAs 
arise in two type of regions: (1) an H° region 
completely opaque to ionizing photons with hv > 
13.6 eV and (2) a mildly ionized region containing 
intermediate-ionization species such as Al 2+ , but 
not high ions such as C +3 or Si +3 . The observed 
column density of the i-th ionization state of the 
element X is given by the relation 

N(X*) = f Xk(K)A k (X)n k (R)dl (1) 

k=l,2 Jk 

where k = 1 and 2 indicates the neutral and ion- 
ized region, respectively; the integrals are carried 
on along the portions of line of sight I through the 
two regions; rife(X) = X^n^X 2 ) i s the local den- 
sity (atoms cm~ 3 ) of X summed over all the pos- 
sible ionization states i; Xk(X l ) = nk(X 1 )/ rifc(X) 
is the ionization fraction of the i-th state; A k (X) 
= rife (X)/rifc (H) is the absolute abundance of X. 

We assume that type-1 and type-2 regions 
have equal abundances in a given DLA system: 
Ai(X) = A 2 (X) = A(X). For the dominant ion- 
ization state in the type-1 region, id, we have 
Xi(X ld ) = 1 and, in particular, £i(H°) = 1. From 
these assumptions and from (1) we can express 
the intrinsic abundance ratio of two elements X 
and Y in terms of the column density ratio of the 
dominant species 



A(X) _ N(X id ) 
A(Y) ~ N(Y^) 



x C(X/Y) , 



where 



C(X/Y) = 



7Vi(H)+ J 2 x 2 (Y^)n 2 (H)rf/ 
JVi(H) + / 2 a:2(X*««)n2(H)<a 



(2) 



(3) 



is, by definition, the ionization correction factor 
and Aq(H) = J m(H)dl the mean hydrogen den- 
sity in the H° region. It is easy to obtain similar 
a expression for the correction term of absolute 
abundances, C(X/H), by replacing Y ld with H° in 
Eq. (3). 

We can also derive the column density ratio be- 
tween two ionization states of a given element. In 
particular we are interested in comparing the ion- 
ization state which is dominant in the neutral re- 
gion, id, with higher ionization states, ih- Since 
xi(X ih ) = and xi(X ld ) = 1 we obtain from (1) 



n 



X th _ JV(X** ) _ J 2 x 2 (X th )n 2 (H)<fl 
XV = N(X* d ) ~ JVi(H) + J 2 x 2 (X^)n 2 {U)dl ' 

(4) 

We can derive simpler expressions by introduc- 
ing the average ionization fractions along the line 
of sight, x 2 {X l ) = j 2 x 2 (X*)n 2 (tt)dl/N 2 (R). With 
this definition we obtain 



and 



x 2 (Y^) + 
x 2 {X id ) + N 1 /N 2 

x 2 (H°) +Nx/N 2 
x 2 {X^) + Ni/N 2 

X^ = x 2 (X^) 
{ X^' x 2 {X i d) + N 1 /N 2 



C(X/Y) = 



C(X/H) 



(5) 



(6) 



(7) 



where Nx/N 2 = Ni (H) /N 2 (H) is the fraction of 
the line-of-sight total hydrogen column densities 
in the two regions. The mean ionization frac- 
tions a; 2 (X l ) can be estimated by modeling the 
intensity and spectrum of the ionizing radiation 
field in which the clouds are embedded. The 
parameter N\/N 2 plays a central role in assess- 
ing the importance of ionization effects. When 
Ni/N 2 ^> 1 the ionization corrections are neg- 
ligible since C(X/Y) ~ 1 and C(X/H) ~ 1. In 
this case the ionization ratio TZ(X th /X ld ) tends to 
be very low no matter which are the conditions 
in the ionized envelope. On the other hand, if 
N 2 /N\ » 1, the ionization corrections and the 
ionization ratio are representative of the ionized 
envelope and not of the neutral region. This lat- 
ter case is appropriate when we consider that all 
species of low and intermediate ionization state 
originate in a single, partially ionized layer. 
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Table 1 

Al+ and Al 2+ column densities in DLA systems. 



QSO 


Z abB 


log JV H i 




log AT(A1 2 +) 


log AT(A1+) 


Rof 


log JV(A1+) 


d a 


!°S X', 2 + + ? 

° N(Al+) 


0000-263 


3.3901 


21.41 


± 


0. 


08 


12.54 ± 0.10 b 


> 13.15 


2 


13.66 


± 


0. 


.17 


-1.12 ± 0.20 


0100+130 


2.3090 


21.40 


± 





05 


12.72 ± 0.03 


— 


1 


14.17 


± 





25 


-1.45 + 0.25 


0149+33 


2.1400 


20.50 


± 


0. 


10 


12.56 ± 0.04 


12.94 ± 0.10 


1 


13.14 


± 





18 


-0.38 ± 0.11 


0201+365 


2.4620 


20.38 


± 


0. 


04 


13.61 ± 0.01 


- 


5 


14.14 


± 





.24 


-0.53 ± 0.17 




1 7688 


20.00 


± 





18 


1 -3 on _|_ n n7 




2 


13.46 


± 





20 


— u.zo zc u.z± 


0216+080 


2.2931 


20.45 


± 


0, 


16 


13.74 ± 0.02 


> 13.88 


2 


14.03 


± 





22 


-0.29 ± 0.22 


0307-4945 


4.466 


20.67 


± 


0. 


09 




13.36 ± 0.06 


7 


13.28 


± 


0. 


.17 




0458-02 


2.0400 


21.65 


± 


0. 


09 


13.34 ± 0.02 


> 13.72 


1 












0528-2505 


2.1410 


20.70 


± 


0. 


08 


12.77 ± 0.08 


> 13.46 


2 


13.84 


± 





19 


-1.07 ± 0.20 


0528-2505 


2.8110 


21.20 


± 


0. 


10 


> 14.07 


> 14.20 


2 












0841+129 


2.3745 


20.95 


± 


0. 


09 


12.60 ± 0.10 c 


> 13.26 


1 


13.82 


± 





.18 


-1.21 ± 0.21 


0841+129 


2.4764 


20.78 


± 


0. 


10 


12.63 ± 0.04 


> 13.19 d 


1 












0951-0450 


3.8570 


20.60 


± 


0, 


10 




13.30 ± 0.02 


1 


13.22 


± 


0. 


.17 




1104-1805 


1.6616 


20.85 


± 





01 


13.06 ± 0.02 


>13.43 d ' e 


3 


13.96 


± 





21 


-0.90 ± 0.21 


1215+33 


1.9990 


20.95 


± 


0, 


07 


12.78 ± 0.02 


> 13.36 


1 


13.61 


± 





17 


-0.84 ± 0.17 


1331+170 


1.7764 


21.18 


± 





04 


12.97 ± 0.01 f 


> 13.57 


1 


13.86 


± 





19 


-0.88 ± 0.19 


1346-0322 


3.7360 


20.70 


± 


0. 


10 




12.55 ± 0.03 


1 


12.52 


± 


0. 


.17 




1425+603 


2.8268 


20.30 


± 


0, 


04 




> 13.55 


2 












1795+75 


2.6250 


20.80 


± 


0. 


10 


13.62 ± 0.04 f 




1 


14.11 


± 





24 


-0.49 ± 0.24 


1946+765 


1.7382 










12.70 ± 0.04 




2 


13.33 


± 





17 


-0.63 ± 0.18 


1946+765 


2.8443 


20.27 


± 





06 


< 12.09 


12.26 ± 0.03 


2,6 


12.19 


± 





.17 


< -0.17 


2206-199A 


1.9200 


20.65 


± 


0, 


10 


13.88 ± 0.01 f 




4 


14.39 


± 


0. 


.28 


-0.51 ± 0.28 


2206-199A 


2.0760 


20.43 


± 


0. 


10 




12.18 ± 0.01 


4 


12.32 


± 





.17 




2230+025 


1.8642 


20.85 


± 


0, 


08 


13.60 ± 0.01 f 


> 14.02 


1 


14.24 


± 





.26 


-0.64 ± 0.26 


2231-0015 


2.0662 


20.56 


± 


0, 


10 


13.14 ± 0.03 c 




1 


13.83 


± 





.19 


-0.69 ± 0.19 


2237-0607 


4.0803 


20.52 


± 





11 




12.85 ± 0.02 


2 


12.84 


± 





17 




2348-147 


2.2790 


20.56 


± 





08 


11.95 ± 0.08 s 


12.65 ± 0.01 


1 


12.79 


± 





.17 


-0.70 ± 0.08 


2359-0216 


2.0950 


20.70 


± 


0, 


10 


12.88 ± 0.05 8 


> 13.66 d 


1 


13.99 


± 





.23 


-1.11 ± 0.24 


2359-0216 


2.1540 


20.30 


± 


0, 


10 


12.96 ± 0.02 


13.17 ± 0.02 


1 


12.88 


± 





20 


-0.20 ± 0.03 



a Value estimated from 7V(Si~*~) and from the linear relation shown in Fig. 1. The error bar is obtained from error 
propagation of the 7V(Si + ) error and of the linear regression error. The latter is estimated by taking into account 
the dispersion of the linear regression, a, and the errors <J m and o q given in the third row of Table 2. 

b Value determined from the UVES spectrum presented by Molaro ct al. (2000a). Only the Al 2+ line at 186.2 
nm has been used. The line at 185.4 nm is affected by telluric contamination. 

d The published value is treated as a lower limit since the line is saturated in a substantial part of the profile. 

c Mcan value of the lines at A 185.4 and 186.2 nm. 

c The profile fitting result given by the authors is not reliable owing to the strong saturation of the line; we adopt 
as a lower limit the result of the optical depth analysis. 

f The Al 2 ~*~ absorption profile shows one or more features which are not seen, or are much weaker, in the absorption 
profiles of low ions; the equivalent width of such features is significantly smaller than that of the bulk of the 
absorption seen also in low ions. 

s The column density refers only to the Al 2 "^ absorption profile in the radial velocity range where the bulk of low 
ions is observed; a high-velocity Al 2+ absorption component is also present, but its column density is not reported 
here. 

References.— (1) Prochaska & Wolfe (1999); (2) Lu ct al. (1996); (3) Lopez ct al. (1999); (4) Prochaska &c 
Wolfe (1997); (5) Prochaska & Wolfe (1996); (6) Lu ct al. (1995); (7) Dessauges-Zavadsky et al. (2001) 
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Table 2 

Linear regression analysis of logarithmic column densities in DLA systems. 



Ions 


n 


r 


a 


m ± a m 


q±a q 


Si+ vs. Fe+ 


30 


0.96 


0.17 


0.99±0.06 


0.64±0.82 


A1+ vs. Fe+ 


8 


0.91 


0.20 


1.20±0.22 


-3.73±3.10 


A1+ vs. Si+ 


9 


0.94 


0.16 


1.05±0.15 


-2.20±2.11 


Al 2 + vs. Fe+ 


20 


0.51 


0.37 


0.46±0.18 


6.25±2.70 


Al 2 + vs. Si+ 


17 


0.74 


0.35 


0.83±0.19 


0.45±2.94 



Table 3 

Ionization correction terms for abundance ratios in DLA systems' 1 



log AT(H°) 


20.2 


20.4 


20.6 


20.8 


21.0 


21.2 


21.4 


21.6 


21.8 


C[N/0] 


+0.002 


+0.001 


+0.001 


+0.001 


+0.000 


+0.000 


+0.000 


+0.000 


+0.000 


C[N/Si] 


+0.149 


+0.100 


+0.066 


+0.043 


+0.028 


+0.018 


+0.011 


+0.007 


+0.004 


C[N/S] 


+0.227 


+0.157 


+0.105 


+0.070 


+0.045 


+0.029 


+0.019 


+0.012 


+0.007 


C[N/Fe] 


+0.027 


+0.017 


+0.011 


+0.007 


+0.004 


+0.003 


+0.002 


+0.001 


+0.001 


C[N/Zn] b 


+0.642 


+0.498 


+0.372 


+0.269 


+0.188 


+0.128 


+0.085 


+0.055 


+0.036 


C[0/Mg] 


+0.152 


+0.102 


+0.068 


+0.044 


+0.028 


+0.018 


+0.011 


+0.007 


+0.005 


C[0/Si] 


+0.147 


+0.099 


+0.065 


+0.042 


+0.027 


+0.017 


+0.011 


+0.007 


+0.004 


C[0/S] 


+0.225 


+0.155 


+0.105 


+0.069 


+0.045 


+0.029 


+0.018 


+0.012 


+0.007 


C[0/Fe] 


+0.025 


+0.016 


+0.010 


+0.006 


+0.004 


+0.003 


+0.002 


+0.001 


+0.001 


C[0/Zn] b 


+0.640 


+0.496 


+0.371 


+0.268 


+0.187 


+0.127 


+0.085 


+0.055 


+0.036 


C[Mg/Si] 


-0.005 


-0.004 


-0.002 


-0.002 


-0.001 


-0.001 


-0.000 


-0.000 


-0.000 


C[Mg/S] 


+0.073 


+0.053 


+0.037 


+0.025 


+0.017 


+0.011 


+0.007 


+0.004 


+0.003 


C[Mg/Fo] 


-0.127 


-0.086 


-0.057 


-0.037 


-0.024 


-0.015 


-0.010 


-0.006 


-0.004 


C[Mg/Zn] b 


+0.488 


+0.394 


+0.304 


+0.224 


+0.159 


+0.109 


+0.073 


+0.048 


+0.031 


C[A1/Si] 


-0.531 


-0.431 


-0.334 


-0.249 


-0.178 


-0.123 


-0.082 


-0.054 


-0.035 


C[A1/Fe] 


-0.653 


-0.514 


-0.389 


-0.284 


-0.201 


-0.137 


-0.092 


-0.060 


-0.039 


C[A1/Zn] b 


-0.038 


-0.033 


-0.028 


-0.023 


-0.017 


-0.013 


-0.009 


-0.006 


-0.004 


C[Si/Fo] 


-0.122 


-0.083 


-0.055 


-0.036 


-0.023 


-0.015 


-0.009 


-0.006 


-0.004 


C[Si/Zn] b 


+0.493 


+0.397 


+0.306 


+0.226 


+0.160 


+0.110 


+0.074 


+0.048 


+0.031 


C[P/Si] 


-0.150 


-0.111 


-0.079 


-0.054 


-0.036 


-0.024 


-0.015 


-0.010 


-0.006 


C[P/Fe] 


-0.272 


-0.194 


-0.134 


-0.090 


-0.059 


-0.039 


-0.025 


-0.016 


-0.010 


C[P/Zn] b 


+0.343 


+0.286 


+0.227 


+0.172 


+0.124 


+0.086 


+0.058 


+0.038 


+0.025 


C[S/Si] 


-0.078 


-0.057 


-0.039 


-0.027 


-0.018 


-0.011 


-0.007 


-0.005 


-0.003 


C[S/Fe] 


-0.200 


-0.139 


-0.094 


-0.063 


-0.041 


-0.026 


-0.017 


-0.011 


-0.007 


C[S/Zn] b 


+0.415 


+0.341 


+0.267 


+0.199 


+0.143 


+0.098 


+0.066 


+0.044 


+0.028 


C[Ti/Si] 


+0.083 


+0.057 


+0.038 


+0.025 


+0.016 


+0.010 


+0.007 


+0.004 


+0.003 


C[Ti/Fe] 


-0.039 


-0.026 


-0.017 


-0.011 


-0.007 


-0.004 


-0.003 


-0.002 


-0.001 


C[Ti/Zn] b 


+0.576 


+0.454 


+0.344 


+0.251 


+0.176 


+0.120 


+0.08 


+0.052 


+0.034 


C[Cr/Si] 


+0.134 


+0.091 


+0.060 


+0.039 


+0.025 


+0.016 


+0.010 


+0.006 


+0.004 


C[Cr/Fc] 


+0.012 


+0.008 


+0.005 


+0.003 


+0.002 


+0.001 


+0.001 


+0.001 


+0.000 


C[Cr/Zn] b 


+0.627 


+0.488 


+0.366 


+0.265 


+0.185 


+0.126 


+0.084 


+0.055 


+0.035 


C[Mn/Si] 


+0.086 


+0.059 


+0.039 


+0.026 


+0.017 


+0.011 


+0.007 


+0.004 


+0.003 


C[Mn/Fe] 


-0.036 


-0.024 


-0.016 


-0.01 


-0.006 


-0.004 


-0.003 


-0.002 


-0.001 


C[Mn/Zn] b 


+0.579 


+0.456 


+0.346 


+0.252 


+0.177 


+0.121 


+0.080 


+0.053 


+0.034 


C[Ni/Si] 


+0.042 


+0.029 


+0.020 


+0.013 


+0.009 


+0.005 


+0.004 


+0.002 


+0.001 


C[Ni/Fe] 


-0.080 


-0.054 


-0.035 


-0.023 


-0.015 


-0.009 


-0.006 


-0.004 


-0.002 


C[Ni/Zn] b 


+0.535 


+0.427 


+0.326 


+0.239 


+0.169 


+0.115 


+0.077 


+0.050 


+0.033 


C[Zn/Fo] 


-0.615 


-0.480 


-0.361 


-0.262 


-0.183 


-0.125 


-0.083 


-0.054 


-0.035 



a Results predicted by model S2 at log U — — 2.2 (see Section 3.1). All values are given in logarithm. 
b Systematic errors probably present owing to uncertainty of zinc atomic parameters 
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log N(Si + ) 

Fig. 1. — Comparison of Al + and Si + column den- 
sities in DLAs. Continuous line: linear regression 
of the data points. Dashed lines: ±la dispersion 
of the linear regression. 
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log N(H"| 

Fig. 3.— Model predictions of the ft(Al 2 +/Al+) 
ratio versus 7V(H°) in DLAs in the case of a soft, 
stellar ionizing source (model S2 discussed in Sec- 
tion 3.1.1). Solid curves: log7e(Al 2+ /Al+) pre- 
dicted at \ogU = —2.6 (bottom curve), —2.2 
(thick line), and —1.7 (top curve). Diamonds: ob- 
servational data points as in Fig. 2; dashed-dotted 
line: linear regression through the observed data; 
dotted lines: ±ler dispersion of the linear regres- 
sion. 
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log N(H°) 

Fig. 2. — Ionization ratio 7\L(A1 2+ /A1 + ) versus 
neutral hydrogen column density in DLAs. Data 
points taken from Table 1. 
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Fig. 4.— Model predictions of the ft(Al 2 +/Al+) 
ratio versus _/V(H°) in DLAs in the case of 
a hard, QSO-dominated ionizing source (model 
HI discussed in Section 3.1.2). The continu- 
ous curves have been obtained at constant value 
of the photoionization parameter, U. Curves 
from bottom to top correspond to logU = 
-4.8,-4.2,-3.6,-3.0,-2.4,-1.8 and -1.2, re- 
spectively. Filled and empty symbols represent 
solutions for which 72.(Si +3 /Si + ) < —0.5 dex and 
ft(Si +3 /Si+) > -0.5 dex, respectively. Dashed- 
dotted line: linear regression through the observed 
data points. Dotted lines: ±lcr dispersion of the 
linear regression trough the observed data. 
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Fig. 5. — Ionization correction terms for [Mg/H] 
measurements in DLAs. Continuous curves: pre- 
dictions for the S2 model discussed in Section 3.1.1 
calculated at constant values of the ionization pa- 
rameter U. Thick curve: \ogU = —2.2. Thin 
curves: \ogU = —1.7 (bottom) and \ogU = —2.2 
(top). Empty circles: predictions for the HI model 
discussed in Section 3.1.2 calculated at logU = 
—4.2 (left symbol) and log U = —4.8 (right sym- 
bol). 
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Fig. 8. — Ionization correction terms for [P/H] 
measurements. Legend as in Fig. 5. 
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Fig. 9. — Ionization correction terms for [S/H] 
measurements. Legend as in Fig. 5. 



Fig. 6. — Ionization correction terms for [Al/H] 
measurements. Legend as in Fig. 5. 
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Fig. 7. — Ionization correction terms for [Si/H] 
measurements. Legend as in Fig. 5. 
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Fig. 10. — Ionization correction terms for [Ar/H] 
measurements. Legend as in Fig. 5. 
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Fig. 11. — Ionization correction terms for [Ti/H] 
measurements. Legend as in Fig. 5. 
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Fig. 12. — Ionization correction terms for [Cr/H] 
measurements. Legend as in Fig. 5. 
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Fig. 14. — Ionization correction terms for [Fe/H] 
measurements. Legend as in Fig. 5. 
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Fig. 15. — Ionization correction terms for [Ni/H] 
measurements. Legend as in Fig. 5. 
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Fig. 13. — Ionization correction terms for [Mn/H] 
measurements. Legend as in Fig. 5. 
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Fig. 16. — Ionization correction terms for [Zn/H] 
measurements. Legend as in Fig. 5. 
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